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We study the role of quark-quark correlation in the baryon structure and, in particular, 
the hyperon non-leptonic weak decay, which is sensitive to the correlation between quarks in 
the spin-0 channel. We rigorously solve non-relativistic three-body problem for SU(3) ground 
state baryons to take into account the quark-pair correlation explicitly. With the suitable 
attraction in the spin-0 channel, resulting static baryon properties as well as the parity 
conserving weak decay amplitudes agree with the experimental values. Special emphasis is 
placed also on the effect of the SU(6) spin- flavor symmetry breaking on the baryon structure. 
Although the SU(6) breaking effects on the local behavior of the quark wave functions are 
considerable due to the spin-0 attraction, the calculated magnetic moments are almost the 
same as the naive SU(6) expectations. 



§1. Introduction 

Properties of light baryons have been extensively studied by various models based 
on the constituent quark picture, in which the constituent quarks are assumed to be 
identified with quasi-particles of non-perturbative QCD vacuum. Their results for 
static hadron properties are consistent with experiments including applications for 
the two nucleon systems, although this model involves several adjustable parameters. 
Despite the success of this approach, it is not clear whether or not the constituent 
quark model correctly describes the quark distributions in the baryons, and further 
provides the non-leptonic weak hyperon decay with AI = 1/2 rule. 

The AI = 1/2 rule implies dominance of the AI = 1/2 transitions and strong 
suppression of AI = 3/2 process in the non-leptonic hyperon decay This em- 
pirical rule does not originate from the theory of the weak interaction itself, and 
therefore one needs an explanation of some dynamical origin. As we will show in 
section 2, the non-leptonic weak decay of hyperons is reasonably described by the 
soft-pion theorem. It leads to the baryon pole approximation which can reproduce 

*' E-mail: hiyama@post.kek.jp 
**' E-mail: ksuzuki@la.numazu-ct.ac.jp 
***■' E-mail: toki@rcnp.osaka-u.ac.jp 
****' E-mail: kami2scp@mbox.nc.kyushu-u.ac.jp 



typeset using pTpTEX.sty <ver.0.8> 



2 



E. Hiyama, K. Suzuki, H. Toki and M. Kamimura 



relative magnitudes of various hyperon decays and AI = 1/2 rule. Within this ap- 
proximation, the weak decay takes place as the two-quark transition process where 
a US-pair in the initial hyperon with their total spin changes to a spin-0 ud-pair 
in the final state baryon. However, if one calculates its matrix clement using the 
constituent quark model, the absolute value of the amplitude is about a half of the 
experimental data at most *^ . We emphasize here, because of the heavy T^-boson 
mass, the weak matrix elements are quite sensitive to the short range quark-quark 
correlations. Consequently, the failure of the constituent quark model to describe 
the non-leptonic weak transition may suggest the lack of the quark correlation in 
the s = channel. Several works ^^'^^ also suggest the importance of the short- 
range quark correlation on the AI = 1/2 decay. In particular, it is known that the 
parity-conserving mr~^ decay is free from factorization and penguin contribu- 

tions, and simply dominated by the two-quark transition process in the baryon pole 
term . Therefore, this decay mode gives a crucial constraint on the strength of the 
quark-quark correlation. 

On the other hand, it was recently pointed out from both theoretical^^' and 
phenomenological points of view, there exists a strong correlation between quarks 
in the s = channel, corresponding to the correlation between the quark-antiquark 
spin-0 pair which forms the pion as the highly collective state. The simple constituent 
quark model has never incorporated such a quark-quark correlation properly. 

In this paper we try to clarify the role of the quark-quark correlation in the 
baryon structure and hyperon non-leptonic weak decays by using the constituent 
quark model. We assume the short range spin-dependent correlations between quarks 
together with the confinement force, and calculate the baryon masses and other static 
properties. In order to deal with the spin-dependent correlation in the correct and 
systematic way, we must rigorously solve the three-body problem. For this purpose, 
we adopt the Gaussian expansion method (GEM) for few-body systems, which has 
been developed by the two of the present authors (E.H. and M.K.) and their collab- 
orators " (see Ref . 9 for a review) . We assume the isospin symmetry between u 
and d quarks, and solve the three-body problem without further approximations or 
assumptions. Our calculations do not rely on the SU(3) fiavor symmetry, and hence 
the strange quark is distinguished from light u, d quarks. To our knowledge, this 
work is a first attempt to study the non-leptonic hyperon decay and other hadron 
properties consistently in the framework of the constituent quark model by taking 
into account the quark-quark correlation. 

So far the constituent quark model with the confining force and the perturba- 
tive gluon exchange provides reasonable results for masses and radii, but fails in 
reproducing the hyperon weak decay After introducing the suitable quark-pair 
correlation, one does not know a priori whether or not the constituent quark model 
can explain both standard properties and at the same time non-leptonic weak decay 
matrix elements. One of the main purposes of this paper is to make this issue clear 



*' It may bo possible to reproduce the experimental values by artificially taking a smaller baryon 
radius. Such a procedure cannot be justified, because the baryon radii should be also chosen to be 
consistent with experiments. 
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by solving the non-relativistic three-quark problem explicitly. 

We also focus on the SU(6) breaking effects on baryon properties. Introduction 
of the spin-dependent correlation naturally spoils the SU(6) spin- flavor symmetry 
which is known to work well for e.g. the baryon magnetic moments. Hence, we 
shall calculate the magnetic moments to estimate the SU(6) breaking effects. We 
discuss how the SU(6) breaking affects the baryon properties, and point out that 
this symmetry is still useful for the static baryon properties, although local behavior 
of the quark wave function considerably departs from the SU(6) symmetric limit. 

This paper is organized as follows. In section 2, we introduce the effective weak 
interaction which includes the renormalization group improved QCD corrections. 
We show the relevant formulae within the pole approximation, and define a set 
of matrix elements. In section 3 we construct the non-rclativistic potential model 
which incorporates the confinement force and the quark correlation. Our numerical 
procedure to solve the three-body problem is described in detail here. In section 4 we 
calculate the baryon masses, radii, and the non-lcptonic weak transition amplitudes. 
Comparison with experiments will be made there. We also evaluate the magnetic 
moments of the SU(3) baryons in order to clarify the SU(6) breaking effects in section 
5. Final section is devoted to the summary and discussions. 

§2. Calculation of weak matrix elements 

We write the low energy effective weak interaction Hamiltonian density ; 

^, , , Gp sm9 cos6 sr-^ , ^ , ^ , , ^ 

nw{x)= ^ ^Q(//2)0i(x)+h.c. , (2-1) 

where 

01 = [u7^(l - J5)s][d7''{l - 75)ti] 

02 = [(^7^(1 - I5)s][u7''{l - ^5)u] 

q=u,s,d 

04= E [97,.(l-75)s][dr(l-75)g] 

q=u,s,d 

05 = [d>(l-75>] E [^"7^(1 + 75)9] 

q=u,s,d 

06= E [97^(1 -75)s][dr(l + 75)9]. 

q=u,s,d 

This effective weak Hamiltonian can be obtained by integrating out over ly-boson 
degrees of freedom in the Standard Model. Current-current operators Oi , O2 give 
dominant contributions in our case, and O3 ~ Oq provide so called penguin contri- 
butions. The coefficients Ci^fx"^) of the Hamiltonian get QCD radiative corrections 
calculated by the renormalization group technique, and thus depend on the scale. 
We take the scale = IGeV^ as a typical scale of light hadrons, and use values of 
Ci given in ref. 15). 
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Our task here is to evaluate the matrix element {Bf 7r"|'Hvy(0)|-Bj) for the 
strangeness changing process Bi — > Bf + tt"". We note that the factorization as 
well as the penguin contributions are too small to reproduce the AI = 1/2 ampli- 
tudes. QCD radiative corrections to the effective weak Hamiltonian tend to increase 
the AI = 1/2 amplitudes, but effects are not enough. 

Analysis based on the chiral dynamics of low energy QCD is suitable to deal 
with the strongly interacting pion-nucleon system. With the help of the soft pion 
theorem, one can rewrite the transition matrix element as 

{Bf{pf) 7r"(g)|W^(0) \B,{p,)) = j d^x e-''^-{-q^ + ml) {Bf |T {7r«(x), Wh^(O)} \B,) 

I I d^xAl{x),rLw{^) \B, 



1% 



/ d''xe-"^^{Bf |r{^^(x),Wiy(0)} \B,). (2-2) 

The first term is called the commutator term which gives the parity violating S-wave 
amplitudes. The second term expresses the baryon pole contribution by inserting 
the intermediate baryon states between j4^(x) and TLwi^), ^"^^ contributes to both 
parity conserving and violating amplitudes, as illustrated in Fig. 1. Here, the initial 



Bj ^ Bn Bf Bj Bn ^ Bf 



~7l 



71 



□= 




Fig. 1. Baryon pole contributions to the non-leptonic hyperon decay. Thick and thin solid lines 
denote the baryon and quark, respectively. The pion is depicted by the dashed line. The weak 
transition between baryon states is shown by the filled box, which expresses the two quark weak 
transition process us — > ud shown in the right hand side. 



hyperon Bi changes to an intermediate state baryon Bn by the weak interaction and 
then Bn emits the pion to produce the final state i?y(and vice versa). In this work, 
we concentrate on the parity conserving P-wave amplitudes to investigate the quark 
correlation inside the ground state baryons*^ For example, the parity conserving 
amplitude for yl'^ ^ n + vr*^ is given by 



B{A^ 
Mn + Ma 

U 



n + TT 



0\ 



G 



1 



''Ma - Mn 



{n\npc\A) + {n\npc\r 



1 



Mn - Me 



GttO 
AS 



(2-3) 



where {n\T-Lpc\A) and {n\T-Lpc\^^) are the matrix elements between appropriate 
baryon states with Tipc being the parity conserving part of the weak Hamiltonian 

*' The parity violating amplitudes are not adequate for the study of the quark correlation 
selectively, since the contributions from the commutator, penguin, factorization, and the baryon 
pole are all comparable. 
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defined in the Appendix. G^^, denotes the axial vector couphng constant which 
gives a probabihty for the pion emission B ^ B' + -jt°'. Here, we consider only the 
ground state baryon octet as the intermediate states. Formulae for other hyperon 
decays are found in the Appendix. 

The axial vector coupling constants are rather well-known quantities from ex- 
periments. We adopt the values of G™^/ obtained by the SU(3) parameterization 
for the existing experimental data, since the axial vector coupling seems to be insen- 
sitive to the quark correlation due to its one-body operator structure. Therefore, we 
are now in the position to determine the matrix elements of the weak Hamiltonian, 
(nlHpcl^) and {p\'Hpc\S'^)*'' ■ We recall that the quark models such as Isgur-Karl 
Harmonic Oscillator model or MIT bag model give much smaller values for 
these matrix elements than the experimental data^^'^^. It is instructive to rewrite 
the weak Hamiltonian TLpc in the non-rclativistic limit in the coordinate space as 

sing cos g jvj^ 
npc = ^ (ciOi +C2O2 ), (2-4) 

Of ^ = -02^« = 4 4(1 - au ■ ds) S^'^ {rus) a^as , (2-5) 

where Oj, a| are annihilation and creation operators of quarks with flavor i. Presence 
of the spin-projection operator {I — au ■ ffs) in Oi^, tells us that only the 
(ns) pair with their total spin being can contribute to the weak decay process, 
namely, the weak transition is generated by the two body process between spin-0 
quark pairs; (us)^ (ud)^. The isospin of the initial {us)^ pair is 1/2, and the final 
(tid)" has the isospin-0 due to the antisymmetrization. Thus, this process guarantees 
the AI = 1/2 dominance in the non-leptonic hyperon decays as pointed out long 
ago"*"^. Now it is clear that this decay amplitude is very sensitive to the correlation 
of the spin-0 quark pair in the baryons'^\ The standard constituent quark model 
never incorporates such a correlation properly. However, in fact, some fundamental 
studies on non-perturbative QCD^^'^^ suggest that there exists the strong attractive 
correlation for the quark-quark pair with s = 0. These considerations naturally lead 
us to study the quark structure of baryons by taking into account the attractive 
correlation which could enhance the weak decay amplitudes. 

§3. Constituent quark model for baryons with spin-dependent 

correlation 



Our purpose here is to construct the quark model to deal with the quark pair 

correlations and thus account for the non-leptonic weak decay. Although there are 
several efforts from the lattice QCD simulation and phenomenological analysis, our 
knowledge of the interaction between light quarks is still far from complete under- 
standing. As we have discussed in the introduction, the instanton liquid model of the 

non-perturbative QCD vacuum provides a spin-dependent correlation between the 
quark-pair. Due to the finite spatial size of the instantons, typically 0.35fm , such a 

*^ In this paper, we restrict ourselves to study only A, E hyperon decay. Results including other 
hyperons will be published in subsequent publication. 
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spin-spin interaction should not be point-like but has a finite range. The constituent 
quarks are also assumed to have their internal structure. Hence, we introduce the 
Gaussian shape spin-dependent force which acts on only the quark-quark pair with 
s = 0. We neglect possible flavor dependence of the potential. Other spin-dependent 
pieces like the spin-orbit and tensor interactions are neglected for simplicity, because 
we shall clarify the effects of the spin-spin interaction on the non-leptonic weak decay 
amplitudes. 

On the other hand, we take the two-body Harmonic Oscillator potential as the 
confinement force, since the analytical solutions of the three-body system are well 
known In order to check accuracy of our numerical calculations, we can compare 
analytical results with ours, when we turn off the spin-dependent interaction. Choice 
of the HO potential is advantageous for us to develop our numerical procedure in 
this paper, but it can be easily improved in a more realistic way. 

Finally, we phenomenologically introduce the effective Hamiltonian which in- 
cludes the confinement force Vc and the spin-dependent part Vg as 

H = Y,^-Tg + Vc + Vs + Vo, (3-1) 



^C = E^^(^^-^^)'' (3-2) 

i<j 

Vs = Y -^^exp [- (xi - Xjf (spin=0 pair) (3-3) 

and V5 = for spin=l pair. Here, mj, Xj and pj are the mass, coordinate and 
momentum of the z-th constituent quark, and K,Css,P are the model parameters 
which are taken to be common for all the baryons concerned. Tg is the cm. kinetic 
energy and Vq is the constant parameter which contributes to the over all shift of 
the resulting spectrum and is chosen to adjust the energy of the lowest state to the 
nucleon mass. The quark masses are taken to be m„ = = 330 MeV and = 500 
MeV. 

Using this Hamiltonian, we shall solve non-relativistic three-body problem rig- 
orously. We assume only the isospin symmetry between up and down quarks. The 
quark wave functions are constructed by the antisymmetrization without invoking 
any further approximations or assumptions. We note that the SU(6) spin-flavor 
symmetry is broken within our formalism because of the spin-dependent correlation. 
Namely, the spatial part of the total wave function is dependent on the spins and 
isospins of the three quarks and is expanded in terms of a number of basis functions 
so as to describe the spin-dependent short-range correlations. 

Since A^, A and f2~ are composed of three quarks having the same isospin and 
the other baryons A,U,E*S and H* are not, we introduce two different types of 
total wave functions for the two cases. 
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3.1. Wave functions ofN,A and 

According to the Gaussian expansion method " , we consider three rearrange- 
ment Jacobian coordinates of Fig. 2 and refer them as channels c = 1 to 3; here, 
Tfc = Xj — Xj and = x^ — (xj + Xj)/2 for the cyclic permutations of {i,j, k). We 



1 

Q 



1 



1 




2 '1 3 
c=1 



Fig. 2. Rearrangement Jacobian coordinates for N, A, and Q 

first construct three-body basis functions for the spin, isospin and spatial part of the 
channel c = k with J, M (total spin and its z-component) and T, Tz (total isospin 
and its z-component) as 



c=k) 



[Xi(«) Xi(j)]5Xi(fc) 



JM 



[Vr{i) rir{j)\tr]T{k) 
with the isospin of quarks, r, to be 



TTz 



^ = {s,S,n,l,N,L,I,t} (3-4) 



I {u,d) for AT, /A 
(s) for f2- . 



(3-5) 



Here, t/t-t^ {i) is the isospin function of the i-th quark, t denoting the isospin of each 
quark-pair. Xl^(^) is the spin function of the i-th. quark, s and S denote the intrinsic 
spin of each quark-pair and three quarks, respectively. The numbers n and I {N and 
L) specify respectively the radial and angular-momentum excitations with respect 
to the Jacobian coordinates Vc (Rc), and / and L are coupled to the total orbital 
angular momentum I. Explicit form of the spatial functions 4'ni{^k) and ■0ArL(Rfe) 
will be discussed below. The three-body center-of-mass motion does not appear in 
our framework. 

The totally antisymmetric basis function with the quantum number set ^, ^jmtt^. 

(c) 

is obtained by the equal-weight superposition of '^jmtt^ s, ^^^^ c = 1 — 3, multiplied 
by the color singlet wave function and posed by the Pauli restriction l-\-s-\-2T-\-t-\-l = 
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even : 



^jMTT^,e,= ^(color singlet) 



c=l 



JMTT:,, i 



(3-6) 



The totally antisymmetric property of the basis functions ^jmtt^,^ is explicitly seen 
by interchanging the particle numbers of any pair. 

The total wave function, ^jmtt^, is given as a sum of these basis functions: 



(3-7) 



This form is the most general one of the totally antisymmetric three-quark functions 
of N,A and f2~ (and their spatially excited states). The coefficients ^4^ are to be 
determined by solving the Schrodinger equation 



(H - E) ^JMTT. = 



(3-8) 



with the Rayleigh-Litz variational principle. 



3.2. Wave functions of A, E,E*,^ and ^* 

In this case, in order to construct totally antisymmetric wavefunctions, we con- 
sider according to the Gaussian expansion method " , the nine rearrangement 
Jacobian coordinates of Fig. 3 (7 = 1 — 3, c = 1 — 3) in which a particle (illustrated 
by a double circle) is s quark and the other two are u,d quarks for A,E,E*, and 
the situation is opposite for S, S*. The channel name 7 indicates the particle num- 
ber of the sole, different quark (double circle), whereas the channel name c denotes 
the particle number which the Jacobian coordinate R points. In order to make the 
coupling scheme of the spins and isospins of the three quarks as visual as possible, 
we always place the sole, different quark as the thirdly (lastly) coupled particle in 
the spin-isospin space (but not always in this order in the coupling scheme of the 
coordinate space). 

The three-body basis function for the spin, isospin and spatial part of the channel 
7, c is given by 



(7, c) 

JMTT„i 



Xi(a)xi(/?)J^Xi(7) 



t>nl{r-i,c)i>NL{^n.c) 



J JM 



rir{a)rir{l3) 



^ = {s,S,n,l,N,L,I,t} (3-9) 



where a, /?, 7 are given by the cyclic permutations. The isospins r and r-y are defined 
as 



i {u,d) foryl,^,^* 
(s) forH,H*, 



is) 



, , f^l' ^'5(3-10) 

^ {u,d) for , ^ ^ 



Using these basis functions, we can construct two types of totally antisymmetric 
three-quark basis functions: 
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G- 



2 '11 3 
(Y,c)=(1,1) 



Rl2A''l2 ""is A Ri 







2 3 
(Y.c)=(3,3) 





2^3 
(Y.c)=(3,1) 




2 h^ 3 




2 3 
(Y.c)=(3,2) 



A-Type 



B-Type 



Fig. 3. Rearrangement Jacobian coordinates for A, S, S* , ^ , and 



(i) A-type: equal-weight superposition of the three basis functions with 7 = c = 
1 — 3 (cf. the left-most column of Fig. 3), multiplied by ^(color singlet) and posed 
by the Pauli restriction l-\-s-\-2T-\-t-\-l = even for the quark pair having the same 
isospin (strangeness): 



^jMrT.,«= ^(color singlet) )_J ^'J^'tt^^^ ■ 

7=c=l 



Fi(7,c) 



(3-11) 



(ii) B-type: equal-weight superposition of the six basis functions with 7 7^ c = 
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1 — 3 (cf. the middle and the right columns of Fig. 3), multiplied by <P(color singlet): 

= ^(color singlet) j^E'^JMTT.,^ ■ (3-12) 

7=1 

The total wave function is given by a sum of these basis functions: 

This is the most general form of the totally antisymmetric three-quark wave functions 
of A, and S* (and their spatially excited states). The A-type basis functions 

are useful for describing the correlations between the quarks having the same isospin 
(strangeness) whereas the B-type ones are effective for the correlations between the 
quarks which have different isospins (cf. Fig. 3). The coefficients ^4^ and B^i are 
determined by the variational principle. 

It is to be noted that, when calculating the energy E and the coefficients and 
5^/, use of the channel 7 = 1 alone in Eq. (3-12) and (3-13) is sufficient since the 
strong interaction does not mix the configurations having different 7. But, use of 
the full wave function with 7 = 1 — 3 (replacing the so-obtained A^ by A^/\/?t and 
B^, by S^//\/3) is necessary in the calculation of the weak decay matrix elements. 

3.3. Spatial basis functions 

Our wave function allows complicated admixture of spin-isospin states depending 
on the quark-pair correlations. For the spatial basis functions ^nZm(r) and t/jNLMi^), 
we have to employ any functional form which satisfies the following requirements: i) 
The functions should be very suited for describing both the short-range correlations 
and the long-range tail behavior, ii) Energy matrix elements should be calculated 
analytically and easily between the basis functions of the different arrangement set 
of Jacobian coordinates, iii) Non-linear parameters of the basis functions can be 
searched quickly. 

To the authors' opinion, the most suitable are the Gaussian basis functions with 
range parameters chosen to lie in a geometrical progression ^) " : 

4>nim{r) = Nni r'e-('-/'^")V,^(r), 
V'iVLM(R) = Nnl e-W^^)'yLM(R), (3-14) 

where N^i and N^l are normalization constants and 

l^n — fmin(^ (jT' = 1 f'-max)) 

RN = RminA''-^ {N=l^Nr^^). (3-15) 

Successful applications of the Gaussian expansion method with the use of the above 
Gaussian basis functions are seen in Refs. 9)- 14) for various three- and four-body 
systems. 
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3.4. Matrix elements of the Weak Hamiltonian 

With the use of the above wave functions, the weak decay matrix elements, 



qNR 



E+ 



n 



QNR 



n 



QNR 



Firstly, we consider the operation of t' 



Aj, are calculated in the following manner: 
le operator Eq. (2-5) on the hyperon wave 
function. Since s quark changes to u quark and d quark to u quark due to the 
weak interaction (Fig. 1), the operation on the isospin part of the wave function 
results in 



Qd'^i (1 - • o-s) <5(x„ -Xs) Qutts I r?i i(q;)??i 1 (/3)r/oo(7) ) 
= (1 - <Ta • (Ty) (5(xa -x^) I r/i_i(a)r?i i(/3)r/i 1(7) ) 
+ {l-(Tf3- a^) S{yip - x^) I r?i i(a)r/i_i(/?)r7i 1 (7) ) 

The same operation on the isospin part of the A wave function gives 



(3-16) 



44 (1 



1 

fa • 



r]ii{a)r]i_i{P) - ?7i_i (a)r/i 1 r/oo(7) ) 



(5(xc 



x^ 



??i_i(a)??i_i(/?)??i 1(7) ) 
22 22 



+ - '^13 ■ 0-7) <^(x/3 - x^) I ??i i(q:)??i 1 1(7) ) • (3-17) 

•t / ^ 2 2 2 2 2 2 

Overlap between the proton wave function and the so-operated wave function, 
multiplied by Gi7sin^cos0(ci — C2)/\/2, gives the amplitude {p\Hpc\ E~^), and sim- 
ilarly for {n\Tlpc\A). From a simple relation between Clebsch-Gordan coefficients, 
we have (n |7ipc| U<^) = {p |7ipc| /V2. 

§4. Mass spectrum and structure of the baryons 

Wc shall fix the model parameters K.,(3 and Cgs/'nif^ so as to reproduce the 
baryon masses and the charge radius. The experimentally measured proton charge 
radius includes contributions from both valence quark core part and its meson clouds. 
It is reasonable to subtract the vector meson dominance contribution from the data 
of the proton electric charge radius (0.86fm)^ to obtain the valence quark core ra- 
dius (rp)core- From this analysis, we extract (r^) 

core ~ (0.6fm)^. By searching the 
parameters within a reasonable range, we obtain the parameters K = 0.007GeV^, 
P = 0.55 fm and Css/ml ~ l.lOGeV which give m{A) - m{N) ~ 293 MeV and 
(^p)corc = (0.60fm)^. After this determination, there are no more adjustable pa- 
rameters in our calculation. With these parameters, we obtain the neutron charge 
radius square (r^) = — 0.05fm^, which is also consistent with the experimental value 
(r^) = — 0.12fm^, after taking into account contributions from the meson cloud. 

As for the angular momentum space and the Gaussian range parameters of the 
basis functions, we examined that the following case is good enough. Contribution 
of the orbital angular momenta other than Z = L = 7 = is found to be negligible 
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as long as the calculation is made of the ground states of baryons without non- 
central forces as in the present study. As an example, spin s and S and Gaussian 
range parameters employed for the case of nucleon are listed in Table I only for 
I = L = I = 0. The same range parameters are used commonly for any set of spins 
and isospins of other baryons and are suitable enough to obtain good convergence 
of the eigenenergies. 

Table I. Three-body angular-momentum space and Gaussian range parameters (in fm) for the 
ground state of nucleon (J = 1/2+, T = 1/2). 





L 


I 


s 


S ^max 


^min 


r N 

' max " max 


^min 


-^max 














1/2 8 


0.1 


3.0 8 


0.1 


3.0 











1 


1/2 8 


0.1 


3.0 8 


0.1 


3.0 



Wc show first in Fig. 4 SU(3) baryon mass spectrum. The agreement with the 
data encourages us to proceed with our approach. It could be possible to adjust the 
model parameters or to elaborate the form of the potential to get much improved 
spectrum, but the present results are enough for our purpose here. 

[MeV] 
1800 



1600 



1400 



1200 



1000 
800 

Fig. 4. SU(3) baryon mass spectrum: Calculations are shown by the solid lines, and experiments 
by the dashed ones. 

To clarify the effects of the attractive correlation in the spin-0 quark pair on the 
nucleon structure, we introduce the quark-pair correlation function. 

pW(r) = (^,^,^J<5(x, - X,- - r)p(-)(zi)|lP-i^,i^J, (4-1) 

where r stands for the distance between a quark pair and the P^*^ (ij) the projection 
operator of the quark-pair spin s; in other words, the quantity p^'^^r) is the proba- 
bility density to find a spin-s quark pair at the distance r. The density correlation 
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function at origin p^^\0) in the s = case essentially fixes the value of the weak 
decay matrix element in Eq. (2.3). p^^\r), which is independent of the angle f, is 
illustrated in Fig. 5. In Fig. 5 we find a large deviation between the s = and s = 1 
density distributions. p^^\0) in the s = case is about three times as large as that 
of the s = 1 case. Similar tendency is seen in other baryons. This provides a huge 
enhancement for the AI =1/2 weak decay amplitude. 

0.4| ■ < ■ < ■ < ■ < ■ 1 

\ 

\ s=t=1 

\ 

s=t=0 

7^ \ 
E ■ 

H — \ 

£ 0.2- \ 




r (fm) 

Fig. 5. Quark-pair density distribution p''''(r) with the spin-dependent correlation. Results for 
s = t = 1 and s = t = pairs are depicted by solid and dashed curves, respectively. 

For the nucleon, we calculate the r.m.s distance f^^^ and R^^^ with respect to 
the Jacobian coordinates r and R by f^*) = [/ r'^ p^^\r)dr]^/^ and similarly for R^^\ 
We get f(*=°) = 0.92 fm and ^(^=°) = 0.97 fm, while f(*=i) = 1.10 fm and ^(^=^) = 
0.81 fm. Apparently, the s = correlation considerably contacts f and extends R 
compared with the s = 1 case, but the r.m.s matter radius of the total system is 
almost 0.58 fm both for the s = and 1 cases. It seems that the correlation is not 
so strong to form the so called 'diquark'-clustering in the nucleon^). 



§5. Numerical results for weak decay amplitudes 



With the parameters fixed in the 
previous section we calculate the ma- 
trix elements of the weak Hamiltonian 
shown in Table II. In the left column 



Table II. Matrix elements of the operator, 
Of ^ (in 10"^GeV^) with and without the 
spin-dependent attraction Vs. 

with Va without Va 
(n|Of -0.960 -0.399 

(p|Of^|r+) 2.760 0.977 
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we show the matrix elements with the 
quark correlation and those without the 
correlations in the right column. In the 
absence of the correlation V5 = 0, £1 Tel- 

tio {p\Oi^\I!^) / {n\Oi^\A) = —2.45 shows a perfect agreement with the SU(6) 
expectation — \/6 ~ —2.4494- • •. This agreement ensures the validity of our numeri- 
cal calculations (Not only the ratio but also the absolute value have been examined). 
In the realistic case with the spin-dependent force, one can observe the substantial 
enhancement of the matrix elements and the SU(6) breaking effects. 

The non-leptonic weak transition parity-conserving amplitudes are tabulated in 
Table III. We show the pole contributions only in the second column and additional 

Table III. P-wave non-leptonic weak transition amplitude (in lO"'^ unit). In the second and forth 
columns, values in brackets show the results without K -Empirical values are taken from ref. 1). 



Decay 


Pole 


others 


Total 


Exp. 




23.23(13.45) 


2.05 


25.28(15.50) 


26.6 


r+ ^ n7r+ 


40.9(8.21) 


0.00 


40.9(8.21) 


42.2 




-5.29(-3.54) 


-5.02 


-10.31(-8.56) 


-15.8 



factorization and penguin contributions, taken from Ref. 4), arc in the third column. 
Then, we show the total decay amplitudes in the forth column to be compared with 
the experiments. It is worth noting that the i7+ — > nTr"*" parity conserving decay 
process is completely dominated by the baryon pole diagrams without any factor- 
ization or penguin contributions. This fact indicates that the Z""*" mr'^ P-wave 
amplitude is the most appropriate observable to probe the quark-quark correlation 
in the baryons. We find a good agreement for E — Nir decays, while the pole 
contribution of the A Ntt is not enough. In general, our calculations reasonably 
reproduce the magnitudes of the pionic hypcron decay amplitudes. However, there 
is a large cancellation between the first and the second terms in Eq. (2-3), which 
makes the resulting A — > Ntt amplitude small. Hence, the pole contribution to 
A Ntt process strongly depends on the values of G^^, G^n well as the weak 
matrix elements. Slight variation of G]^ ^, ~ 10% yields about 40% modification 
of the A — ^ Ntt pole contribution, whereas the decay amplitudes for i7 — > A*" -|- vr" 
are essentially unchanged. The A decay amplitudes are also rather sensitive to the 
detailed shape of the potential. 

We have obtained the agreement when we take the Gaussian size parameter 
P = 0.5 ~ 0.6 fm for the spin-dependent potential. If we choose smaller values of 
the size parameter /3, e.g. 0.2 fm, the weak transition matrix element is enhanced 
and overestimates the experimental data by order of magnitudes. In this work, 
we have retained the non-relativistic kinematics so far. Since the mass of the con- 
stituent quark is comparable with the kinetic energy, we should consider relativistic 
corrections to both the wave function and the calculation of the weak decay matrix 
elements. Inclusion of such relativistic corrections also changes the value of size 
parameter /3. The study of such a dependence is in progress. 

In this paper, we concentrate on the calculation of the parity-conserving P- 
wave amplitudes, since they are particularly sensitive to the correlation strength. 
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On the other hand, for the parity violating S-wave amphtudcs, contributions from 
commutator, factorization, penguins and the negative parity baryon pole terms are 
comparable. This fact suggests that the S-wave amplitudes are not suitable to quan- 
tify the quark-quark correlation. Nevertheless, it is possible to calculate the S-wave 
amplitude with our wave functions. Results for the commutator contributions ampli- 
tudes are also consistent with the data after including the penguin, factorization and 
negative baryon pole contributions. Such a study will be discussed in forthcoming 
paper. 

§6. SU(6) symmetry breaking effects on the magnetic moment 

In the previous sections wc have evaluated the pionic weak transition amplitudes 
as well as the mass spectrum. Due to the strong spin-dependent attraction our wave 
function enhances the weak decay amplitudes and violates the naive SU(6) spin- 
flavor symmetry. As we have shown in Table II, the ratio of the matrix elements of 
the weak Hamiltonian {p\Oi / {n\Oi ^\A) shows a clear evidence for the SU(6) 

breaking. This ratio becomes —2.45 without the correlation in a agreement with the 
SU(6) result — \/6, while —3.02 with the quark correlation. Therefore, we estimate 
the size of the SU(6) breaking effects to be about 20%, which is significant. 

On the other hand, it is historically known that the light baryon magnetic mo- 
ments are well reproduced in the naive quark model by virtue of the SU(6) spin-flavor 
symmetry. Hence, it is also interesting and important to estimate effects of the SU(6) 
breaking by calculating the baryon magnetic moments. 

Neglecting contributions from the quark orbital angular momentum, the opera- 
tor for the baryon magnetic moment is simply given by 



where Qi is the electric charge operator of the i-flavor quark, and i runs over i = 
1 ~ 3. The magnetic moment of the baryon B is obtained by directly calculating 
the matrix element 



in terms of our wave functions. This quantity certainly depends on the internal 
spin-flavor structure of the baryon through mj, Qi and af. 

Our results are shown in Table IV. In the first column, we show the results 
with the spin-dependent force, and those without the correlations in the second 
column which are the same as the naive SU(6) results. It is manifest that the results 
are almost unchanged even after introducing the spin-dependent correlations. The 
differences are of order of a few % in any cases. 

These results tell us that the global spin-flavor structure of quark wave functions 
(integrated over volume) is quite insensitive to the existence of the quark-quark cor- 
relations. As a result, the success of the SU(6) symmetry for bulk baryon properties 




(6-1) 



= {B\flmag\B) 



(6-2) 
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Table IV. Baryon magnetic moments 



with Va without Vs Exp. 



Hp 

Hn 

HA 
He+ 
He- 

Mho 
Hb- 



2.78 2.84 2.792847 

-1.83 -1.90 -1.913042 

-0.602 -0.613 -0.613 ± 0.004 

2.69 2.73 2.458 ± 0.010 

-1.05 -1.06 -1.160 ±0.025 

0.817 0.836 

-1.410 -1.449 -1.250 ±0.014 

-0.507 -0.502 -0.6507 ± 0.0025 

-1.84 -1.84 -2.02 ±0.05 



seems to be maintained. On the other hand, the introduction of the quark correla- 
tion modifies the local structure of the quark wave function substantially shown in 
Fig. 5, and thus enhances the weak decay matrix elements. 

The violation of the SU(6) spin-flavor symmetry in the local quark distributions 
can be found in other experiments, namely, the deep inelastic lepton scattering off 
the nucleon. Measured ratio of the momentum distribution functions d{x)/u{x) at 
X ~ 1 tends to in contrast to the SU(6) value 1/2. This is one of the examples that 
demonstrates the breakdown of the local SU(6) spin-flavor symmetry. We also note 
that such a flavor dependence of the quark distribution functions is also explained 
by considering the quark-quark correlation in the s = channel . 



In conclusion, we have studied the role of the spin-dependent quark-quark cor- 
relation in the baryon structure. We have pointed out that the non-leptonic weak 
transition of the hyperon is an unique quantity to investigate the quark correla- 
tion in the spin-0 channel. In particular, E'^ mr'^ P-wave decay is free from 
the factorization or penguin contributions as shown in Table III, and thus serves 
a severe constraint on the quantitative understanding of the spin-0 correlation. In 
order to demonstrate its importance for the weak transition, we have introduced the 
non-relativistic constituent quark model with the spin-dependent attraction. We 
have solved three-body problem rigorously using Gaussian expansion method ' . 
Such a spin-dependent interaction may originate from the non-perturbative QCD 
dynamics. Results for static baryon properties reasonably agree with the empirical 
values. 

We have developed the procedure to calculate the matrix elements of the weak 
Hamiltonian with a number of the rearrangement channels as discussed in section 3. 
Same technique can be applied to the strangeness — 2-system, S hyperons, and such 
a study is in progress. Calculated transition amplitudes without the spin-dependent 
quark correlation agree with the SU(6) predictions perfectly. This result guaran- 
tees the accuracy of our numerical calculation. Resulting parity-conserving weak 
transition amplitudes are consistent with the experiments, when the quark correla- 
tion is turned on. At the present, quantitative description of the non-leptonic weak 
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hyperon decay is still difficult in any models of the baryons^^'^^^'^^^ Naive chiral 
perturbation theory cannot explain the parity-conserving and violating amplitudes 
simultaneously, and convergence of the chiral expansion seems to be worse . Here, 
we present a possible improvement of this long standing problem in the framework of 
the constituent quark model, but still have several things to work out. For example, 
the A Ntt transition amplitudes are underestimated due to the large cancella- 
tion in the pole formula, although they are strongly dependent on the choice of the 
parameters. 

Introduction of the quark correlation, which is strong enough to explain the non- 
leptonic weak decay, modifies the valence quark structure of the baryons considerably. 
The distance between quarks in the spin-0 pair becomes shorter by 20%. Although 
the modification is not so large to induce the diquark-clustcring of the nucleon, such 
a tendency is consistent with the several phcnomcnological studies 

Although the spin-dependent correlation certainly violates the SU(6) spin-flavor 
symmetry for the wave functions, we have found that the baryon magnetic moments 
are almost unchanged compared with the prediction of the naive SU(6) model, as 
demonstrated in Table IV. This is because the magnetic moments are quite insen- 
sitive to the variation of the wave function at short distances r, i? ~ 0. There is 
a historical argument that the SU(6) symmetry assumption in the quark model is 
indispensable to keep the impressive agreement of the magnetic moment with the 
data. However, it is now evident that, even after including the SU(6) spin-flavor 
breaking effects on the wave function, one can obtain quite reasonable values for the 
baryon magnetic moments. 

Recently, strong evidence for a new five-quark baryon state 0^(1540) (known as 
a pentaquark) has been reported by several groups ' . In Ref . 27) , it is suggested 
that the baryon is a bound state of four quarks and an antiquark, containing two 
highly correlated ud pair. The quark-quark correlation potential obtained in the 
present work to explain the weak decay of hyperons as well as the baryon mass 
spectrum would be useful in the study of the structure of the five-quark baryon. 
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Appendix 

We give some formulae to calculate the weak decay amplitudes. First we sum- 
marize the expressions for the baryon pole contributions to the parity conserving A 
and E weak decay. 



B{A n7r°) 



Mn + Ma 



''''Ma-Mn 



{n\Hpc\A) 



+ (n\Hpc\E' 



1 



Mn - Me 



AE 



(A-1) 



B{2J^ —^ mr^) = 
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where Tipc is the parity conserving part of the weak Hamiltonian; 

^, Gpsin^cos^ - i-.y \,- „ \ /T \/- /; M 
ripc = (ci [{d'y^u){u-f^'s) + (d7/,75'«)(^t7 75*)] 

+C2 [(fi7^n)((i7^s) + (il7M75^i)(d7^75s)]) , (A-4) 

where we use only the current-current operators Oi , O2 which give dominant con- 
tributions here. We omit the penguin contributions, since we are interested in the 
absolute magnitudes of the weak transition amplitudes. We use the following pa- 
rameters of the weak Hamiltonian, Gp = 1.16639 x 10~^GeV~^, sin0 = 0.219, 
cos e = 0.975, and = 92MeV. 

As for the axial vector coupling constant, we adopt the SU(3) Goldberger- 
Treiman relation following to the standard approach. 



GZ = G;', = ^{f + d)g (A.5) 



V3{Mn + Ma) 
2Ms- 



G'sJrS^ ~ ~^iii:o — oiiJ ^ 9 (A-7) 



where g is the strong nNN coupling constant, and /, d the SU(3) coupling with a 
condition f + d = 1. We use / = 0.38 and d = 0.62 to obtain the numerical results 
presented in this paper. As pointed out in the text, the pole contribution A Nit 
is very sensitive to the choice of /, d. 

The matrix elements of the weak Hamiltonian must be evaluated with calculated 
quark wave functions. The isospin symmetry for u and d quarks implies a relation 

n|7ipc|i:°) = ^(p|Wpc|^+) (A-9) 
which exactly holds with our wave functions. 
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